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ABSTRACT 

The Enhanced Tributary Monitoring Program (ETMP) was Initiated in 
1979 and 1980 at 15 key Great Lakes tributary outlet monitoring 
stations to obtain a better data base, primarily by increased 
sampling frequencies; to improve the precision of tributary load 
estimates; and, to establish procedures to optimize future Great 
Lakes tributary sampling. Additional water-quality monitoring was 
undertaken at existing Provincial Water Quality Network (PWQN) 
stations with supplementary fjiding received from Environment Canada 
under the terms of the Canada/Ontaro Agreement. This report 
documents the status of the enhanced monitoring; presents sumnaries 
of the monitoring data; and, reconmends improvements for future 
enhanced monitoring. 

Water-quality measurements and load estimates for the 1980 water 
year (October 1979 to September 1980) are presented for various 
water-quality parameters monitored at the 15 ETMP outlets. The 
precision of annual phosphorus loadings (i.e. expressed as the 
percentage of mean phosphorus load) improved as a result of the 
enhanced monitoring. For total phosphorus, the mean standard error 
for the 15 tributaries improved from 36.2% average for routine 
monthly monitoring to 8.09% average for the more frequent ETMP 
monitoring. Sampling requirements (i.e. frequency and distribution 
in time) for two, principal water-quality parameters (total 
phosphorus and nitrate-nitrogen) were found to be different for each 
monitored tributary. 

Enhanced monitoring has continued to date and its future use will 
continue to allow the ccmputation of precise tributary loads, and 
optimize future resource allocation by further site-specific 
refinements of sampling strategy (i.e. sample-collection frequency 
and distribution in time). 
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1.0 INTRODUCTI 



This report documents the status of the Enhanced Tributary 
Monitoring Program; presents a surmary of the monitoring data; and 
recommends improvements for future enhanced monitoring. 

1.1 Background 

The Ontario Ministry of the Environment (MOE) has been measuring 
water quality at strategic locations throughout Ontario since 1964 
to monitor the effects of point-source waste-water discharges and 
those effects originating frcm various non-point sourcEi land uses. 
As part of the Ministry's contribution to the International Joint 
Commission's programs on the Great Lakes, conputed estimates of 
phosphorus loadings from key Great Lakes tributaries located in 
Ontario are provided to, and reported annually by the Great Lakes 
Water Quality Board. Historically, these phosphorus loads were 
derived from a marginal water-quality data base, as the: Provincial 
Water Quality Monitoring Network (PWQN) lacked the fineincial support 
to obtain frequent water samples during high-flow even1:s. 
Supplementary finding received from Environment Canada under the 
terms of the Canada/Ontario Agreement allowed increased sampling of 
some stations in the existing Network. This sampling was initiated 
in 1979 and the results are discussed in this report. In addition - 
to phosphorus, supplementary water-quality and sediment-quality 
monitoring was undertaken (see Table 2) to enable the canputation of 
other contaminant loads. Sediment-quality measurements are the 
subject of a f uture re pjirt, 

1.2 Objecti ves 

The objectives of the Enhanced Tributary Monitoring Program are to 
obtain a better data base, primarily by increased samp'e 
frequencies; establish procedures to improve the precision of the 
tributary load estimates; and, to establish procedures to optimize 
future Great Lakes tributary sampling. 
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TABLE 1: WATERSHED PHYSICAL CHARACTERISTICS 



I 

ro 
I 



Terminal 
Basin 


Drainage Area 

in Canada 
(103 hectares) 


Tributary T 
Dra 
(10 


ributary 
inage Area 
■^ hectares) 


Percentage of 
Canadian Great 
Lakes Drainage 
Area 


Mean 
1980 


Annual Flow 
Water Year 
cms 


Hydrograph Category 
"natural" or "modified" 
(Section 2.2.2) 


Lake 
Superior 


8.350 


Kaministikwia R. 
Little Pic R. 
Pic River 
Black River 


782 
132 
427 

224 




9.4 
1.6 
5.1 

2.7 






44.2 
16.6 
50.5 
34.0 




modified 

natural 
natural 
natural 




Sub Total 


1,570 




18.8 






145.0 




N/A 


Lake Huron 


8,720 


Mississaugi R. 
Spanish River 
Saugeen River 


930 

1,410 

396 




10.7 

16.2 

4.5 






131.0 

169.0 

63.0 




modified 
modified 
natural 




Sub Total 


2.740 




31.4 






363.0 




H/A 


Lake Erie 


1,600 


Sydenham River 
Thames River 
Grand River 


320 
565 
694 




20.0 
35.3 
43.4 






25.7 
68.2 
77.1 




natural 
natural 
natural 




Sub Total 


1.580 




98.7 






171.0 




N/A 


Lake Ontari 


io 3,090 


We 11 and River* 
Twelve-Mile Or.* 
number River 
Don River 
Trent River 


88.7 
31.6 
1.270 




2.9 

1.0 
41.1 






28.2 

200.0 

8.6 

4.5 

168.0 




modified 

modified 

natural 

modified 

modified 




Sub Total 


1.390 




45.0 






409.0 




N/A 


TOTAL 


21,800 




7.280 




33.4 






1090. 







* Lake Erie diversion via Welland Canal 



TABLE 2: WATER QUALITY PARAMETERS MEASURED AT THE ENHANCED TRIBUTARY MONITORING STATIONS 



Tributary 



Suspended 
Solids 



Total 
Phosphorus 



Filtered Reactive 
Phosphate-P 



Kaministiquia River X 

Little Pic River X 

Pic River X 

Black River X 

Mississaugi River X 

• Spanish River X 

w Saugeen River * 

' Sydenham River X 

Thames River X 

Grand River X 

We 11 and River X 

Twelve Mi 1e Creek X 

Humber River X 

_Dan River X 

Trent River . X 



X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 



X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 



Nitrate- 
Nitrogen 



X - sampled by the contract observers and MOE Regional staff 
* - sampled by the Water Survey of Canada 
NS - not sampled 



X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 



Copper 



X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 



Lead 



Cadmium 



X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 



X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 



Mercury 



PCB PCP 2,4, 5-T 



X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


NS 


NS 


NS 


X 


NS 


NS 


NS 


X 


X 


X 


X 


X 


NS 


NS 


NS 


X 


NS 


NS 


NS 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


NS 


NS 


NS 
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2.0 PROGRAM STUDY APPROACH 

2.1 Tributary Selection 

HistoriexilTy-, total phosphorus loadings have been computed by the 
JIOE^for the IJC Water Quality Board for approximately 60 significant 
Great Lakes tributary streams. Figure 1 illustrates both the 
standard error of phosphorus loadings and their magnitude for these 
streams for the 1973-1978 period. From this group, 15 tributaries 
were selected for inclusion in the Enhanced Tributary Monitoring 
Program. Eleven tributaries with loadings of greater than 300 
kilograms per day were selected on the basis of greatest magnitude 
(Figure 1). Four additional tributaries were selected from a group 
of 12 streams with loading rates in the 100 to 300 kilogram per day 
range. Selection of these four tributaries was governed primarily 
by the desire to attain a more uniform allocation of monitoring 
among the Great Lakes. The tributaries selected are listed in Table 1. 

2.2 Mater Quality 

The water-quality parameters measured at each of the monitoring 
sites are listed in Table 2. Analyses of the water samples were 
conducted primarily at the MOE Central laboratory located in 
Toronto. Some trace elements and nutrient analyses from the Lake 
Superior basin, were conducted at the MOE Thunder Bay laboratory. 

The water-quality data collected in this program are published with 
the data from the remainder of the Provincial Water Quality Network 
(PWQN) in the annual Water Quality Data Series (Volume XV, 1979 and 
Volume XVI, 1980 Water Quality Data. Ontario Lakes and Streams). 

2.2.1 Sample Collection Methods 

Local contract observers were hired to collect water-quality 
samples. The sampling apparatus, collection techniques, 
preservation and storage were consistent with those adopted for the 
PLUARG program (Onn, 1980). Specially prepared sample containers, 
chemical preservation and cold storage were used to inhibit chemical 
reactions and microbial activity. 
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Figure 1 . Variation of the historic phosphorus loading estimates ( 1 973-1 9781 at the six:y Great Lakes 
significant tributary streams. 
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Representative samples were collected by depth-integration 
techniques. Samples were collected by raising and lowering the 
sample container at a constant velocity through one or more sampling 
verticals at each monitoring station. The number of sampling 
verticals varied at each monitoring site from one to three because 
of the seasonal variability in the magnitude of streamflow and the 
vridth of each tributary cross-section. 

2.2.2 Sample Collection Frequency 

To achieve an effective sampling frequency, the 15 tributaries were 
characterized through examination of their hydrographs into two 
basic categories, "natural" and "modified". The "natural" category 
exhibits flow patterns of a high spring runoff period followed by a 
recession to low sunmer flows (see figure 2a). Tributaries in the 
"modified" category exhibit patterns of flow resulting from the 
presence of major dams, diversions and/or urbanization (see figure 
2b, 2c) which are distinctly different from those shown in the 
"natural" category described above. 

Although the basic sample strategy was similar for all sites (I.e. 
frequent sample collection during high flow), two plans for sampling 
frequency were adopted. Sampling frequencies for the monitoring 
sites grouped within the "natural" category were keyed to three 
predetermined flow ranges; high-event flow, event flow and 
baseflow. Sampling frequencies for the monitoring sites within the 
"modified" category varied from once per week to once per day with 
emphasis directed towards periods when high flows were expected to 
recur. Annual sampling frequency for both types of hydrograph 
categories ranged from 60 to 110 samples. 
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Figure 2. Streamflow hydrograph categories. 
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2.3 Water Quantity 

Records of streamflow for this program were generally obtained from 
the Water Resources Branch of Environment Canada. Some additional 
records detailing discharge diversions from the Welland Canal to 
both the Welland River and Twelve Mile Creek were obtained from the 
Western Region Office, St. Lawrence Seaway Authority. 

Where possible, water-quality sampling sites were located coincident 
with established Environment Canada streamflow gauges to facilitate 
the use of mean daily discharge records for the conputatlon of load 
estimates. Additional records of discharge every hour were also 
obtained for major runoff events at these sites. These hourly 
discharges were used to aid in the definition of start and stop 
times for specific events when interpreting the data record. 

In cases where no suitable streamflow gauges exist at sites, 
discharges for the watershed outlets were estimated by using data 
from (one or more) upstream gauges. This was accomplished by 
developing flow estimator equations based on drainage areas. For 
example, on the Sydenham River, water-quality sampling was carried 
out at Wallaceburg, downstream of the confluence of the north branch 
of the Sydenham River (Bear Creek) and the main stream. Gauges are 
located about 20 kilometres upstream of Wallaceburg as follows: 

Station 02GG004 - Bear Creek above Wilkes port - 609 km^ 
Station 02GG007 - Sydenham River near Dresden - 1240 km^ 

Estimates of mean daily discharge at the Wallaceburg site, with a 
total watershed area of 3200 km , were derived with the equation: 

Total flow = flow at 02GG004 x 1.9 + flow at 02GGOO7 x 1.6 

To avoid problems with lag time, no attempt was made to provide 
estimates of discharge e\/ery hour for tributaries requiring the flow 
estimator equation. 
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2.4 Loading Canputation Technique 

The objective of the loading computation is to determine the mean of 
the load population (L) for the required time interval. 

All tributary loading computations included the use of i ratio 
estimator attributed to E.M.L. Beale and advocated for ijse since 
1976 by the International Joint Commission (IJC). This estimator is 
discussed In detail in the IJC Quality Control Handbook for Pilot 
Watershed Studies (International Reference Group on Great Lakes 
Pollution from Land Use Activities, 1975, Section 7.5). The ratio 
estimator Involves the relationship between two related populations, 
in this case flow and load. The flow population Is used to improve 
the estimate of the load by assuming that: 

a) flows are monitored continuously and the mean flow (Q) may 
be determined; 

b) discrete observations are available on the locd population 
(1) v/iere 1. = c.q. for 1. = load, c- = 
concentration, and q. = flow of sample; 

c) the distributions of the flows and concentrations are 
approximately normal. 

Use of the ratio estimator assumes that the populations of flows and 
concentrations are approximately normally distributed; however, the 
annual populations are most often highly skewed. Flows may vary 
over 2 to 3 orders of magnitude while concentrations of suspended 
solids can vary 3 to 4 orders of magnitude. The great^ast variations 

are observed during rainfall-runoff events which, although they 

occupy a relatively small portion of the year on a time basis, 

contribute significantly to the total annual runoff and load. An 
effective means of treating such populations is to divide them into 
relatively homogeneous strata within which the subpopulations more 
closely approach normality. Within strata, the ratio estimator Is 
used to estimate the load and variance. Estimates for individual 
strata are then combined arithmetically to produce a pooled estimate 
of load for the annual period. 
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Analyses were performed with a modified version of a FORTRAN program 
^originally developed by PLUARG investigators. Stratification of the 
record was a subjectiw procedure based on consideration of both the 
annual flow record and the sample concentration data. Flow events 
were identified and the behaviour of concentrations with respect to 
flows noted. The concentrations for most parameters were found to 
hold a generally positive relationship with flow (i.e. increasing 
during flow events in proportion to the magnitude of the event) 
allowing definition of concentration ranges for various ranges of 
flow. Strata boundaries were then defined by assigning flow cut-off 
values. For example, at the mouth of the Grand River, 
suspended-solids concentrations for the springmelt period ranged 
from 230 to 784 mg/L with a mean daily flow of 504 m"^/s. While 
data for all parameters were reviewed, experience showed that 
focusing on the most sensitive parameter, suspended solids, led to a 
stratification scheme suitable for all other parameters without any 
great loss in precision. 

For the tributaries with "natural" flow conditions, a four-strata 
approach as follows was generally used: 

stratum 1 - event flow 

stratum 2 - high-event flow 

stratum 3 - cold water/winter baseflow 

stratum 4 - warm water/sunmer baseflow 

For the tributaries with heavily regulated flows (Kam, Mississaugi, 
Spanish, Welland, and Trent Rivers and Twelve Mile Creek), a 
two-strata approach was generally used. The spring period was 
placed in stratum 2 and the remainder of the year in stratum 1. 
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3.0 DATA INTERPRETATION AND DISCUSSION 

3.1 Water Quality Measurements 

As previously inentioned, water-quality samples collected in the ETM 
program were analyzed at the Ministry of the Environnent's 
laboratories located In Toronto and Thunder Bay. Water-quality 
samples collected from Lake Superior tributaries wer(? analyzed In 
Thunder Bay with the exception of those samples Intended for PCB, 
PCP and 2,4,5-T analyses which were forwarded to the Toronto 
laboratory. All other water-quality samples were delivered to and 
analyzed at the Toronto laboratory. 

The laboratory measurements (i.e. pollutant concentrations) of 
water-quality samples collected during the 1980 water-year (October 
1979 to September 1980) are illustrated for suspended solids; 
nitrate-nitrogen; total phosphorus; reactive phosphate; cadmium; 
copper; lead; mercury; PCB; and PCP respectively. The concentration 
ranges (i.e. minimum to maximum concentration) for each 
water-quality parameter are denoted with a solid vertical line on 
logarithmic scale, for each of the tributary outlet monitoring 
stations. An arrow appears at the lower end of the vertical line 
where one or more of the measured concentrations is less than the 
analytical detection limit reported by the laboratory. The range of 
analytical detection limit concentrations is shown for water-quality 
parameters since detection limits are variable from day-to-day and 
are not consistent between laboratories. The number of samples and, 
where applicable, the percentage of samples exceeding the Provincial 
Water Quality Objective (PWQO) is also shown. The -Mow-weighted 
mean concentration is marked on the vertical line with the symbol 
"X" or where one or more analytical results were fojnd to be less 
than the detection limit, the flow-weighted mean concentration is 
marked on the vertical line with the symbol "<". Wiere sufficient 
data exist, the standard error Is marked on the vertical line with 
horizontal bars, above and below the flow-weighted nean 
concentration. 
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Figure 3a. 1 980 water-year, tributary concentration measurements: suspended solids, nitrate-nitrogen. 
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Figure 3b. 1 980 water-year, tributary concentration measurements: total phosptorus, reactive phosphate. 
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Figure 3c. 1 980 water-year, tributary concentration measurements: cadmium, copper. 
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Figure 3d. 1 980 water-year, tributary concentration measurements: lead, mercury. 
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RgureSe. 1980water-year,tributary concentration measurements: PCB, PCP. 
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As shown by the comparsions of flow-weighted mean concentration and 
concentration range (i.e. minimum to maximum concentra:ion), the 
highest values generally occur in monitored tributaries of lakes 
Erie and Ontario for the variables suspended solids, 
nitrate-nitrogen, total phosphorus, reactive phosphate, copper, lead 
and mercury. Flow-weighted mean concentrations from tiese 
tributaries were used to identify exceedences of the Provincial 
Water Quality Objectives.* The flow-weighted mean concentration Is 
in violation for copper at all tributaries except at the Little Pic 
and Black rivers; and for lead at all tributaries except at the Pic, 
Little Pic, Black and Kaministikwia. The analytical detection limit 
for cadmium is variable (.0001 mg/L to .005 mg/L), both "less than" 
and "greater than" the Provincial Water Quality Objective (.0002 
mg/L) because the analytical methodology was modified during the 
study period. The analytical detection limit for PCB (.02 ug/L) is 
significantly higher than the Provincial Water Quality Objective a 

(.001 ug/L) and PCBs were detected on one or more occasions at all \V 
of the tributaries that were sampled for PCBs (i.e. B" ack, 
Misslssagi and Sydenham were not sampled for PCB). PCP and 2,4,5-T 
are substances with undefined tolerance limits (1978, Ministry of 
Environment) and both compounds were detected on one or more 
occasions at all tributaries that were measured (Blac<, Misslssagi, 
Saugeen, Sydenham, Thames and Trent were not sampled For these 
parameters). 



Current scientific evidence Is insufficient to develop a firm 
Provincial Water Quality Objective for total phosphorus. 
However, in place of a firm objective, the total phosphorus 
concentration of 30 jjig/L is considered as a genei'al guideline 
for rivers and streams. 
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3.2 MuUi -Strata Load Estimates 

Annual load estimates derived from the multi-strata computation 
technique for the variables suspended solids, total phosphorus, 
reactive phosphate, nitrate-nitrogen, cadmium, copper, lead, and 
mercury are provided on Table 3. The table presents mean daily 
loading in kilograms per day. the standard error of the mean and the 
number of samples used to derive the load. 

Concentrations found to be less than the analytical detection limit 
were cormionly measured for the variables cadmium, copper, lead and 
mercury. When this was the case, the detection limit value was used 
for ccmputation and the resulting load estimates are assumed to 
overestimate the true load. Load estimates computed in this manner 
have been flagged with a "less than" (L) remark. 

Figures 4a and 4b are graphic illustrations of the loading estimates 
for the eight water-quality variables reported in Table 3. The 
variables have been separated Into two distinct groups namely, 
nutrients (figure 4a) and trace elements (figure 4b). 

Suspended-solids analyses provide gross measurements of the amounts 
of particulate matter found in the water. They are of interest, as 
suspended solids serve as a transporting agent for a host of_ other 
parameters including some important forms of the nutrients 
(phosphorus and nitrogen) and many of the trace elanents, all of 
v*iich tend to exhibit similar patterns of increasing concentration 
with increasing flow. Examination of the distribution of 
suspended-solids loading estimates (figure 4a) shows 81% of the 
enhanced tributary load is delivered by four tributaries - the 
Sydenham, Thames. Grand and Pic rivers. Together, the Sydenham, 
Thames and Grand rivers, all tributary to Lake Erie, deliver 45% of 
the total enhanced tributary, suspended-solids load. While solids 
enter watercourses from many sources, including sewage treatment 
plants and urban runoff, the intensive agricultural activity in 
southwestern Ontario is likely the cause of the high component in 
the Lake Erie basin. The Pic River in the Lake Superior basin 
provides the greatest individual solids component with 25% of the 
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TABLE 3: 1980 WATER YEAR. TRIBUTARY LOAD ESTIMATES 



Key: 
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Number of Samples 
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Standard Error of the Mean 
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N 
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N 
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62 
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(10.000) 
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(18.3) 
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(3.06) 


55 


4.00 
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11.1 
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56 


0.190L 
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47 


0.083L 

(.003L) 
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Ri«r 
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21 
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(11.6) 


35 
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49 
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(0.07) 


31 
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22 


56. 5L 
(O.OO) 


19 


120L 
(5.79L) 


21 


339L 
(0.00) 


20 


0.300L 

(0.019L) 


- 
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Spanish River 


32 


99,000 
(18.500) 


24 


2390 

(402) 


33 


338 
(32.6) 


27 


37.7 

(5.64) 


23 


72. 7L 
(0.00) 


22 


147L 
(6.9a) 


24 


426L 

(9.64L) 


22 


0.534L 

(0.119L) 


29 


0.311L 

(.016L) 


Saugeen River 


21 


212,000 
(12,300) 


33 


5880 
(422) 


96 
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(19.4) 


33 


90.2 
(8.22) 
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(0.22L) 
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r 1 


1 TA 
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!■> Ql 

[0.45L) 


RI 


(1.72L) 


87 


73. 2L 

(0.8a) 


84 
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(0.08L) 
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(0.54L) 
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0.025L 
(0.002L) 
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0.032L 
(.006L) 
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(5870) 
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(23.5) 
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(14.4) 
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(0.70L) 
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(4.53L) 
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(0.00) 
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48 
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(O.OOIL) 
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(0.040L) 


53 


0.049L 
(.009L) 


Trent River 
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138,000 
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44 
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(22.9) 


44 


61.8 
(7.6) 


50 


72. 9L 
(0.36L) 


30 


' 


" 



) 



monitored suspended-solids load for the 15 tributaries. Extensive 
deposits of lacustrine clay, silt and fine sand within the Pic 
watershed probably contribute to the unusually high solids load. 

The variable of major interest with respect to loadings to the Great 
Lakes, total phosphorus, exhibits the greatest loading rates in the 
Lake Erie tributaries (figure 4a). Together, the Sydenham, Thames 
and Grand rivers deliver 52% of the total phosphorus load for the 
streams monitored in this enhanced program. This large component 
can probably be attributed to the intensive agricultural activity in 
these basins. Ranking next in size. Twelve Mile Creek and the 
Welland River deliver an additional 15% of the load. A major 
component of the flow here is comprised of Lake Erie waters diverted 
to these watercourses by way of the Welland Canal (Table 1). The 
higher loading rates of the Welland River and Twelve Mile Creek are 
a reflection of the large water diversions from Lake Erie as well as 
higher phosphorus concentrations in the Welland River (figure 3b). 
The Pic River of Lake Superior ranks next, delivering approximately 
7% of the total phosphorus load monitored in the enhanced 
tributaries. High phosphorus loads here are attributed to high 
sediment loads probably reflecting natural soil phosphorus losses. 

Reactive phosphate is primarily the dissolved form of inorganic 
phosphorus and is generally considered to be readily available for 
aquatic plant growth. The distribution of load estimates for this 
variable (figure 4a) is very similar to that shown by total 
phosphorus. The three streams tributary to Lake Erie are of major 
importance delivering 59% of the reactive phosphate load monitored 
from the 15 tributaries, followed by the Welland River and Twelve 
Mile Creek with an additional 15%. The remaining 26% is more or 
less uniformly distributed among the other 10 tributaries. 

Nitrate-nitrogen is the end product of the stabilization of organic 
^nitrogen. Elevated levels of nitrate are usually found in polluted 
waters that have undergone some degree of self-purification and also 
occur in watercourses draining fertilized agricultural areas. 
Examination of the distribution of filtered nitrate loading 
estimates (figure 4a) shows 90% of the enhanced tributary load is 
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Figure 4a. 1980 water-year, tributary loading estimates: suspended solids, 
nitrate-nitrogen, total phosphorus, reactive phosphate. 
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Figure 4b. 1980 water-year, tributary loading estimates: cadmium, copper, 
lead, mercury. 
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delivered by six tributaries - the Saugeen. Sydenham, Thames, Grand 
and Welland rivers and Tvgelve Mile Creek. The intensive 
agricultural activity in southwestern Ontario is like'y the source. 
By contrast, the four Lake Superior tributaries deliver only 1% of 
the nitrate-nitrogen load. 

The distribution of loading estimates for the four trace elanent 
variables, cadmium, copper, lead and mercury, show th.it 60 to 70% of 
the load is consistently delivered by four tributaries;, the 
Mississagi, Spanish and Trent rivers and Twelve Mile Creek 
(figure 4b). These four tributaries have an order of magnitude 
higher flows than the remaining 11 tributaries and consequently the 
loads reflect the high volume of water transported past their 
respective outlets. 
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3.3 Improved Precision of Load Estimates 

A major goal of the Enhanced Tributary Monitoring Program (ETMP) is 
to improve the precision of the load estimates through additional 
water-quality sampling to allow the application of a more precise, 
stratified, loading computation technique. Historically, total 
phosphorus loading estimates have been derived from a limited data 
base of 12 to 14 water-quality samples per year. While the historic 
computations have included the use of a ratio estimator, the number 
of water-quality samples has generally been insufficient to allow 
stratification to reduce the bias from the skewed distribution of 
the data. 

In order to evaluate the improvement in precision resulting from the 
current multi-strata load estimates, single stratum load estimates 
were computed for each tributary by selecting 14 samples from the 
ETMP data base (i.e. two samples per month for two high-flow months 
and one sample per month for the remainder of the year) and the 
results compared. This comparison was undertaken for total 
phosphorus, a particulate related parameter, and nitrate-nitrogen, a 
soluble parameter, at each sampling. For certain variables (e.g. 
PCB and cadmium], the loading precision will not improve until 
laboratory detection limits are consistently lower. 

Ccmparison of the current loading estimates computed by multi -strata 
technique with those resulting from a 14-sample, single-stratum 
computation is presented in Table 4. The table presents the mean 
daily loading estimate in kilograms, the standard error expressed as 
a percentage of the mean, and the number of samples used to derive 
the load estimate. For total phosphorus, the mean standard error 
for the 15 tributaries improved from 36.2% for the single stratum, 
14-sample computations to 8.09% for the multi-strata estimates. For 
nitrate-nitrogen, similarly, values improved from 23.2% to 9.55%. 
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TABLE 4: COMPARISON OF MULTI-STRATA AND SINGLE-STRATUM .OAD ESTIMATES: 
TOTAL PHOSPHORUS AND NITRATE-NITROGEN 



V 



Parameters 
Tributaries 


N 


Total Phosphorus 
Multi - N Single- 
Strata Stratum 


N 


Nitrate- 

Multi- 

Strata 


•Nitrogen 
N Single- 
Stratum 


Kam River 
SE {%] 


62 


208 
(8.80%) 


14 


302 

(27.5%) 


51 


336 
(16.2%) 


14 


558 
(35.5%) 


Little Pic R. 
SE 


58 


120 
(9.67%) 


14 


270 

(51.1%) 


21 


ai.8 

(8.04%) 


12 


125 
(41.3%) 


Pic River 
SE 


59 


648 
(7.58%) 


14 


1900 
(62.6%) 


35 


261 
(12.5%) 


14 


435 
(41.8%) 


Bl ack River 
SE 


59 


123 
(8.13%) 


14 


304 
(52.0%) 


34 


318 
(6.10%) 


14 


440 
(31.1%) 


Mississagi R. 
SE 


31 


131 
(9.54%) 


13 


150 
(13.3%) 


25 


1920 

(i:..o%) 


13 


1950 
(18.6%) 


Spanish Ri ver 
SE 


33 


338 
(9.64%) 


13 


354 
(17.0%) 


24 


2390 
(16.8%) 


13 


3180 
(33.3%) 


Saugeen River 
SE 


96 


390 
(4.97%) 


14 


816 
(62.3%) 


33 


5880 
(7.18%) 


14 


6270 
(13.7%) 


Sydenham Ri ver 
SE 


47 


626 
(7.24%) 


14 


844 
(50.6%) 


29 


13.600 
(3,53%) 


12 


14,500 
(20.3%) 


Thames River 

SE 


77 


2190 
(7.99%) 


14 


3200 

(42.8%) 


64 


29,100 
(4,67%) 


14 


27,900 
(12.7%) 


Grand Ri ver 
SE 


117 


1890 
(11.4%) 


14 


1980 
(36.5%) 


56 


16.000 
(4.35%) 


13 


17,200 
(13.4%) 


Well and River 
SE 


105 


468 
(7.69%) 


14 


414 
(20.2%) 


63 


7300 

(9.97%) 


14 


5190 
(16.6%) 


Tv«lve Mile 
SE 


89 


851 
(5.69%) 


14 


1080 
(29.6%) 


43 


6510 
(19.4%) 


14 


5250 
(17.3%) 


Humber River 
SE 


47 


254 
(13.1%) 


14 


128 
(24.9%) 


44 


922 
(9.88%) 


14 


936 
(14.7%) 


Don River 
SE 


61 


256 
(5.63%) 


14 


190 
(14.6%) 


58 


632 
(3.72%) 


14 


760 
(5.45%) 


Trent River 
SE 


66 


535 
(4.28%) 


14 


768 
(37.9%) 


44 


2280 
(7.85%) 


12 


4130 
(31.5%) 


Mean Std. Error 
Std. Deviation 




8.09% 
2.40% 




36.2% 
16.8% 




9.55% 
e.06% 




23.2% 
11.5% 
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3.4 Sample Strategy Analysis 

The required number of samples for various levels of precision can 
be determined theoretically where a prior data base exists and 
strata boundaries have been defined. After satisfactorily 
canpleting the mul ti-strata load estimates, a sampling strategy 
analysis was conducted using an option of the FORTRAN program. 

Assuming sufficient stratification, the number of sanples 
theoretically required to achieve various levels of precision for 
total phosphorus for all 15 tributaries are presented in Figure 5a. 
For example, to achieve +20% precision for total phosphorus sampling 
requirements across the various tributaries range fron 10 (Trent 
River) to 55 (Number River) samples. To Improve the precision to 
+10%, the number of samples theoretically required increases to 44 
and 210 samples. It is clearly evident that a considerable 
expenditure of resources (i.e. additional samples for analysis) 
would be necessary to improve the precision from +20% to +10%. 

Figure 5b presents the same sampling strategy analysis for the 
water-quality parameter nitrate-nitrogen. Here, for +20% precision, 
sampling requirements across the tributaries range from 3 (Sydenhcfn 
River) to 150 (Twelve Mile Creek) samples. To improve the precision 
to +10% the number of samples increases to 12 and 600, respectively, 
again a significant increase. Although not shown, similar results 
could be expected for the remainder of the water-quality variables. 
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Figure 5a. Annual sampling frequency requirements to achieve varying levels of loac ing precision: 
total phosphorus. 
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Figure 5b. Annual sampling frequency requirements to achieve varying levels of loading precision: 
nitrate-nitrogen. 
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4.0 SIMMARY AND CONCLUSIONS 

Supplementary monitoring undertaken at the outlets of fifteen 
Ontario tributaries known to contribute compariti vely high and 
variable phosphorus loads to the Great Lakes has enabled the 
acquisition of a significantly improved water-quality data base. 
This water-quality data base has been used in conjunction with 
streamflow records to compute tributary loading estimates. A 
corresponding improvement in loading estimates is apparent when the 
precision of phosphorus loadings computed with these data are 
compared with phosphorus loads computed from the historic water- 
quality data base. Continuation of the enhanced monitoring will 
provide a high standard of data quality and monitoring uniformity 
(i.e. sampling method, sample collection frequency and water-quality 
parameter coverage) to enable confident future comparisons of 
water-quality measurements and pollutant loading estimates. 

Additional site-specific refinement of sampling strategy (i.e. 
sample collection frequency and water-quality parameter measurement) 
is proposed to optimize resource allocation at the 15 aforementioned 
monitoring stations and to aid in providing for further enhancement 
of future monitoring at two additional Great Lakes tr'butaries, the 
Ausable and Nottawasaga rivers. 
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